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Non-resonant dressing fields can make the transition frequency between two circular Rydberg 
states insensitive to second order variations in the dc electric field. Perturbation theory can be used 
to establish the required dressing field amplitude and frequency. The same perturbative approach 
may be used to understand removal of the first order dependence of the transition frequency on 
electric field about a bias dc electric field [Hyafil et al. Phys. Rev. Lett. 93, 103001 (2004)]. The 
directional alignment of the dressing and dc fields is critical in determining the electric field sensitiv¬ 
ity of the dressed transition frequencies. This sensitivity is significantly larger for circular Rydberg 
states compared to low-angular momentum Rydberg states of Rb. 

PACS numbers: 32.80.Ee, 32.10.Dk, 32.60.-|-i 


I. INTRODUCTION 

Large electric transition dipole moments can exist be¬ 
tween Rydberg states of comparable energies. These 
dipole moments enable strong coupling to electromag¬ 
netic fields at rf/microwave frequencies [1] and enhance 
interactions between Rydberg atoms [2]. Large electric 
transition dipole moments also lead to enhanced dc po¬ 
larizabilities of Rydberg states — the dc Stark shifts of 
low angular-momentum (low-^) Rydberg states scale with 
principal quantum number n like [3]. 

The strong response of Rydberg atom energy levels 
to dc electric fields can be useful, such as for the ac¬ 
celeration and deceleration of the atomic center of mass 
[4], but in other situations this sensitivity is a nuisance 
[5, 6]. For example, as a surface is approached — where 
there are inhomogeneous dc and low-frequency fluctu¬ 
ating electric fields — it may be desirable to maintain 
well-defined Rydberg-Rydberg transition frequencies for 
resonant coupling to surface devices [7]. The large tran¬ 
sition dipole moments of Rydberg atoms enable strong 
resonant coupling, but they also lead to troublesome dc 
Stark shifts, spoiling the resonance condition. 

With the motivation of creating robust qubits, a va¬ 
riety of approaches have been taken in different phys¬ 
ical systems to make resonance frequencies less sensi¬ 
tive to external perturbations. Significant progress has 
been made in superconducting qubits [8] and trapped 
ion qubits [9], to name but two examples. In the case 
of ground electronic state neutral atoms, non-resonant 
“dressing” fields have been shown to reduce the sensitiv¬ 
ity of hyperfine transition frequencies to magnetic field 
fluctuations [10]. 

It is also desirable to reduce the sensitivity of Rydberg- 
Rydberg transition frequencies to low-frequency and dc 
electric fields while maintaining their high sensitivity to 
resonant fields. Hyafil et al. [11] proposed the use of 
non-resonant dressing fields to eliminate the first order 
dependence of the transition energy between two circular 
(jmj=^ = u— 1) Rydberg states on dc field for fluctu¬ 


ations around a non-zero “bias” dc electric field (further 
details are given in Ref. [12]). Bason et al. [13] used non¬ 
resonant dressing fields to increase the sensitivity of Ry¬ 
dberg state energies to dc electric fields, and Sevingli and 
Pohl [14] have explored the influence of multiple dress¬ 
ing fields on Rydberg-Rydberg interactions. Recently, 
Jones et al. [15] experimentally demonstrated a signifi¬ 
cant reduction in the electric field dependence of low-£ 
Rydberg-Rydberg transition frequencies through the ap¬ 
plication of a non-resonant dressing field. 

In the work of Jones et al. [15], the reduction in elec¬ 
tric field sensitivity was discussed in terms of the elimina¬ 
tion of the first-order dependence (dipole) of the dressed 
energy level difference on electric field about a specific 
non-zero dc field: so-called dipole nulling (see Fig. la). 
However, the experimental results showed a stronger sup¬ 
pression of the variation of transition energy with dc elec¬ 
tric field. More specifically, the residual deviations of the 
transition energy were observed to be quartic with dc 
electric field rather than the quadratic behavior expected 
for dipole nulling [11, 12]. In effect, the experimental 
spectra were polarizability nulled (see Fig. la). In con¬ 
trast, the calculations of Hyafil and Mozley et al. [11, 12] 
for circular Rydberg states exhibited only dipole-nulling. 

In the present study, we establish that polarizability 
nulling can be achieved for circular states and compare 
this with the low-^ Rb case [15]. 

Understanding dressed polarizabilities and dipole mo¬ 
ments for circular Rydberg atoms is significant, because 
circular states are both experimentally realizable and 
have special properties such as long radiative lifetimes 
[16]. All atomic species have quantitatively similar cir¬ 
cular Rydberg states due to the low core penetration of 
high-1' states. As shall be shown, the energy degener¬ 
acy of high-1 states gives nulling behavior qualitatively 
distinct from that of low-1 states, particularly when the 
dressing and dc electric fields are not parallel in space. 
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FIG. 1. (a) Polarizability and dipole-nulling contrasted with 
the normal dc electric held dependence of a transition fre¬ 
quency (b) schematic energy level structure near the two 
circular states e and g (m is the parabolic quantum number 
[17]). This diagram is purely qualitative; for example, under 
the conditions studied here, the Stark shifts are much smaller 
than the energy separation between states of consecutive n. 


of the outer (Rydberg) electron to calculate the electric 
dipole couplings. To compactly describe the states of the 
Rydberg electron, we use the approach of Zimmerman et 
al. [18], who consider as a basis set for the dc Stark effect, 
a set of bound states of the Rydberg atom, with energies 
set by the spectroscopy of the zero field Rydberg series. 
With electric fields that produce insignificant ionization 
rates, the continuum states can be neglected. 

In the case of circular Rydberg states — where the 
spin-orbit splitting is insignificant — the basis states are 
labelled as [nfm) where n is the principal quantum num¬ 
ber of the Rydberg electron, £ is its angular momentum, 
and m is the projection of this angular momentum on 
the 2 : axis. 

Due to their limited core penetration, high-£ states 
(such as the circular states) are approximated using hy- 
drogenic wavefunctions. To compute the radial parts of 
the matrix elements between these states, we use Eq.’s 
(63.2) and (63.5) of Ref. [17]. (For these two equations 
to have consistent radial wavefunction phase definition, 
a minus sign should be placed in front of Eq. (63.5).) 

For the \ow-^ states of Rb, spin-orbit coupling can¬ 
not be neglected. Instead we use basis states labeled as 
\ntj'mj) where j and rrij arise from the operator j = s, 
where s refers to the unpaired electron spin. Electric 
dipole couplings between the \ow-i states of Rb are com¬ 
puted in an identical manner to Ref. [18]: the radial Ry¬ 
dberg wavefunctions are integrated in from large r, as¬ 
suming a purely Coulombic potential, with a zero-field 
energy Enij (i.e. Hq \n£jmj) = Enij \n£jmj)) based on 
the known Rb Rydberg state energy levels [19]. Beneath 
a certain low-r the integration is terminated and the 
wavefunction assumed to be zero at lower r. Wavefunc¬ 
tions obtained in this manner can be used to compute the 
radial contribution to {n'£'j'm/l ^\n£jmj). Evaluation 
of the angular contribution to these matrix elements is 
straightforward; we use the formulae of Ref. [18] for 
generalizing for and /Xy. 


II. TECHNIQUES 
A. Dipole matrix elements 

To understand the influence of dressing fields on the dc 
electric field response of Rydberg atoms, we consider the 
zero field Hamiltonian EIq with electric dipole couplings: 

= i7o - A* • Fdc - /X • [(Fac/2) exp(iwdt) + c.c.], (1) 

where fj, is the electric dipole operator and Fdc is the 
dc electric field. The dressing field of angular frequency 
ujd is described using a complex amplitude vector Fac 
which is convenient for the consideration of elliptically 
polarized dressing fields [15]. 

It is assumed that the ionic core of the Rydberg atom is 
in a single state independent of the state of the outer elec¬ 
tron, and that we need only consider the wavefunctions 


B. Eloquet Hamiltonian 

As written, Eq. 1 describes a time-dependent problem. 
However, as the dressing field and thus the Hamiltonian 
is periodic {H{t) = El(t+T), with T = 27r/u;d), Floquet’s 
theorem [20] may be used to write the general solution of 
the time-dependent Schrodinger equation {ihdt |V'(0) = 
E{{t) \ip{t))) using the expansion [21]: 

1x^(1))= (2) 

k 

where the Cfc’s are constant coefficients set by the ini¬ 
tial conditions, the |^fc(f))’s are periodic {\(j)k{t + T)) = 
\4>k{t))), and the Ek's are known as the quasi-energies 
of the dressed system. In what follows, we shall simply 
refer to these as the dressed energies. Since the |(()fe(t))’s 
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are periodic, they may be expanded in a Fourier series: 

\Mt))= E (3) 

9=0,±1,±2,... 

To describe the |^fc(t))’s, we use tensor product notation, 
writing the basis vectors as |nfm) ® |g) = \n£m) 

The periodic Hamiltonian may also be expanded in a 
Fourier series: 

m = E (4) 

p=0,±l,±2,... 

As shown by Shirley [21], the determination of the 
\(j)k{t))’s and EkS that satisfy the time-dependent 
Schrodinger equation amounts to an eigenvalue problem 

HF\Mt)) = Ek\Mt)), (5) 

in which the Floquet Hamiltonian Hp can be written as 

E lEJdq |g) {q\ + E! E{p) 0 |g + p) (g| ■ (6) 

9 9.P 

The Hamiltonian of Eq. 1 allows -Hi? to be written in a 
more specific form: 

Hf =(Ho - pL ■ Fdc) (g) 1 -I- a (g) E |(?) (<?! 

E ® I*?+• 

Q 

( 7 ) 

The matrix elements of Hp in the jnfm) (g \q) basis are 
time-independent. Since we must in principle include all 
possible Fourier components, this matrix is of infinite di¬ 
mension. In practice, Hp can be approximated as a finite 
matrix with a limited number of “sidebands” [22] and di¬ 
agonalized using standard numerical methods [23]. The 
replacement of an infinite matrix representation of Hp 
by a hnite one is a controlled approximation — the size 
of the basis set can be varied and the convergence of the 
eigenvalues and eigenvectors checked. In what follows we 
describe this procedure as a non-perturhative Floquet cal¬ 
culation, to distinguish it from the perturbative approach 
discussed in the next section. 

To characterize the basis sets used for dressed circular 
state calculations, we use three parameters: 6n, Sm, and 
Sq. Writing the specific circular state of interest (either 
e or g) as \nc£ciTic), with £c = nc — £ and TOc = £c, the 
basis set used includes all states \n'£'m') ® \q') with valid 
quantum numbers {£' < n' and jm'j < £') in the range: 
n' = Uc — Sn ,..., ric -f 5n, m' = me — Sm ,... me -I- Sm 
and q' = —Sq ,..., Sq (all ranges are in steps of one). 

When there is only a dc or linearly polarized ac electric 
field present, or they are both present and parallel, the 
quantization axis for the angular momenta is chosen to 
be parallel to the held. In these situations, the Hamilto¬ 
nian matrix separates into diagonal blocks for different 


TO, and the diagonalization effort is signihcantly reduced. 
However in some of the calculations that follow we must 
consider that Fac and Fdc are not parallel, or in general 
that the dressing held is not linearly polarized [15]. In 
these cases, the operator no longer commutes with 
the Hamiltonian, and thus the couplings between states 
of different to must be considered. 

The parameters Sn, Sq and Sm are increased until it 
is observed that the dressed energy — obtained by di¬ 
agonalization in the corresponding basis — is changing 
systematically with basis set size and shows rapid con¬ 
vergence to a hxed value. Unless otherwise stated, we 
have used Sn = 5, Sq = A and Sm = 1. All numerical 
results have been checked for basis set sensitivity. We 
check that an increase in any one of Sn, Sq and Sm from 
the canonical values changes the quoted result by less 
than the precision implied by the number of signiheant 
hgures. 

Similar considerations apply to the basis sets for cal¬ 
culations involving the low-t' momentum states of Rb. A 
formal description is slightly complicated by the non-zero 
quantum defects. For a state of interest \n£jmj), we de¬ 
fine n* using the Rydberg energy: Enij = {—RKh)/n* 
(where -R^b is the Rydberg constant adjusted for the re¬ 
duced mass of e“-Rb''" system). The basis sets used are 
characterized by Sn, Smj, and Sq. They include all states 
\n'£'j'mj') 0 \q') with the zero-field energies of the un¬ 
dressed portion \n'£'j'mj') satisfying 




(n* — SnY 


^ En'Pj' ^ " 


HRb 


SnY 


( 8 ) 


together with mj' ranging over mj — Smj, ...,mj Smj 
and q' ranging over —Sq,...,Sq, both in steps of 1. All 
calculations are done with Sn = 4.5, Smj = 1, and Sq = 
2. The basis size sensitivity of all results are checked in 
a similar manner to the circular states. 

The hyperfine structure of the Rb Rydberg states is 
ignored. Previous experimental work [15] was done with 
the ®^Rb isotope, but all of the results obtained here are 
valid for both abundant isotopes. 


C. Perturbative approach 

Although numerical diagonalization of the Floquet 
Hamiltonian is straightforward, more insight can often 
be obtained by examination of the dominant terms in a 
perturbation expansion of the dressed energies in terms 
of the field strengths. For this purpose we partition Eq. 7 
into a zero-field Hamiltonian: 

Hf,o = Ho 0 1-f 1 (g) E^d? I?) (?! (9) 

q 

and a perturbation term, including both ac and dc fields: 

V =-t E ® {^ac \q + 1) {q\ + Fac* \q - i) (gj} 

9 

-/X • Fdc (g> 1 


( 10 ) 
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such that Hp = Hp^ + V. With basis states of the 
form \n^m) ® |( 7 ), the matrix elements of both V and 
Hp Q do not depend on time and thus time-independent 
perturbation theory may be used to estimate the \4>k{t)) 
and EkS. 

We will initially consider a linearly polarized ac field, 
parallel to the dc field: Fdc = FdcZ, Fac = Fac^, with Fac 
real. Under these conditions non-degenerate Rayleigh- 
Schrodinger perturbation theory may be applied to de¬ 
termine the dressed energies. 

Both the dc and ac Stark shifts are non-vanishing at 
second order in U; to examine the influence of ac fields 
on the dc Stark effect, we must go to fourth order. The 
fourth order shift of the fcth state is (see for example 
[24]): 

A = V (fc| V k) (uj V ju) (uj V jw) (w| V \k) 

w ( 4 °^ - ( 4 °^ - ( 4 °^ - 

\{k\V\u)\^\{k\V\v)\^ , 

u,v ( 4 °) - 

where the summations in u, v and w are over all states 
other than k. We have dropped terms containing the 
matrix elements {k\V \ k) as these couplings vanish in the 
normal zero dc-field Rydberg state spherical basis. The 
terms in this series may be grouped according to their 
dependence on F^c and Fac'- 

= Cfc(4, 0) Fdc^ + Cfc(2, 2, wd) Fdc^ {Fac/2f 

+ C'fc(0,4,a;d) (F,,/2)^ (12) 

where the Ck{a,b,u!dy& have been introduced to repre¬ 
sent terms in groupings by (7^(0, 6, Wd)(Fdc)“(Fac/2)'’. 
For Cfc’s with no dressing frequency dependence (6 = 0), 
LOd is dropped from the argument list. In this notation 
the second order dc Stark shift is /S.Fk = Cfc(2,0)Fdc^, 
whereas the second order ac Stark shift is AEk = 
Cki0,2,iOd)iFac/2)\ 

For polarizability nulling, we are interested in the 
terms that are quadratic in Fdc- 

AUfc [C'fc(2,0)+ ^ Cki2,j,ujd)iFac/2y]. 

7 = 2 , 4 ... 

(13) 

If higher order terms (j > 4) can be neglected, po¬ 
larizability nulling can be obtained when Cfc(2,0) -I- 
Ck{2,2,uJd){Fac/2y = 0. This requires working at a 
dressing frequency such that Ck{2, 0) and Ck{2, 2, ujd) are 
of opposite sign. When there are two states of interest 
(e.g. e and g), the dc differential polarizability is charac¬ 
terized by Ce,g{2^ 0) = Ce{2, 0) — Cg{2, 0), and to obtain 
nulling it is necessary that Ce,g{2, 0) and Ce,g{2, 2, ojd) be 
of opposite signs. 

The dipole matrix elements which contribute to the 
Cfc’s are fixed, whereas the denominators in the pertur¬ 
bation expansion can be varied by tuning the dressing 
frequency. To determine which terms in the perturba¬ 
tion expansion of Eq. 11 are dominant, we examine when 


their denominators show small energy differences (i.e. are 
near-resonant). 


III. RESULTS 


A. Circular Rydberg states with collinear dc and 
dressing fields 


For concreteness — as with Ref.’s [11] and [12] — we 
study the electric field sensitivity of transitions between 
the n = 50 and n = 51 circular Rydberg states (e and 
g hereafter). In the absence of a dressing field, these 
states exhibit a differential dc polarizability characterized 
by Ce,g{2,0) « —25.44Hz/(V/m)^. (Machine readable 
definitions of most calculated quantities and parameters 
in this manuscript are available in Ref. [25].) 

Consider the calculation of Cc{2,2,ujd) using 
Eq. 11 (where c stands for a generic circular state 
\n,£ = n — l,m = n — 1)). Due to the energy degener¬ 
acy of high-t' states with the same n, all three energy 
denominators in the first summation term of Eq. 11 
are resonant when |w) = |n -I- 1, n, n — 1) 0 [g = —1) 
|u) = |n-|-1, n — 1, n — 1) (8) I? = —1) and |u) = \w)- 
For this term the values of 4°^ — Ew\ 4°^ — 4°\ 
and 4°^ — 4°^ are the same (= S) because states of 
the same n but different £ are energy degenerate — a 
situation that does not occur for the lower-^ states of 
non-hydrogenic systems. 

This single triply resonant term with a 1/6^ depen¬ 
dence will dominate over other terms in the series (1/6^, 
1/(5) as the resonance condition is approached. Note that 
the second set of summations in Eq. 11 cannot contribute 
any 1/5^ terms to C'c(2,2, w^). With uJd approximately 
resonant with the n —>■ n-l-1 transition frequency u)c,g we 
have (in atomic units): 


Cc{2, 2 , uje^g F 5) 


\{uh\p-z |wb)P KwbiMz \ch)? , 


where the matrix elements are between the “bare” atomic 
states i.e. Iwb) = |n -I- 1, n, n — 1), etc. In Fig. 2 we show 
Ce,gl2, 2, Ud) as a function of frequency, calculated using 
this simplified form (with two terms, each correspond¬ 
ing to either e or g) and the general series summation of 
Eq. 11. As can be seen, these single terms dominate the 
behaviour near resonance, greatly simplifying our under¬ 
standing of this problem. 

The triply resonant condition found for circular states 
(Eq. 14) is exceptional — for low-£ states there are only 
doubly resonant terms in the perturbation expansion and 
these do not change sign as we go through a resonance. 
Normally, the dressed contribution to the polarizability 
has the same sign on both sides of a resonance - we 
are “stuck” with whatever the resonance gives us. As 
discussed later this is illustrated by the \ow-£ Rydberg 
states of Rb [15]. In contrast, we can make a circular 
Rydberg state either more or less polarizable by sitting on 
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FIG. 2. Coefficient for the nulling of the differential polar¬ 
izability between the two circular states n = 50, £ = 49, 
m = 49 and n = 51, I = 50, m = 50. Note the discon¬ 
tinuity in the vertical axis. The positive (7e,g(2,2, at 
dressing frequencies between the two resonances (n = 51 —>■ 
52 and n = 50 —>■ 51) allows the differential dc polariz¬ 
ability to be nulled by a dressing field with Fac satisfying 
C'e,g(2, 0) -I- <76,9(2, 2,uid){Fac/2Y = 0. “Simplified” refers to 
the use of Eq. 14. 



Fdc (V/m) 


FIG. 3. Resonance frequencies of transitions from the g-state 
(solid lines) and e-state (two dashed lines) as function of dc 
electric field, over the approximate frequency range where 
Ce,g(2, 2, ojd) is positive (see Fig. 2). The two dashed lines cor¬ 
respond to e n = 52, ni = 0,1 (both m = 51) transitions. 
Two specific frequencies (dot-dash) discussed in the text are 
labelled by their detunings 5/27r from the n = 50 —>■ n = 51 
transition frequency. As discussed in Section IIIB, the role 
of orthogonal field fluctuations are investigated at the points 
marked with O. It was not possible to find a dipole null at 
the dc field and dressing frequency marked with A. 


either side of the n —>■ n + 1 transition. Considering that 
(7c (2, 0) is always negative, the absolute polarizability of 
a circular Rydberg state with principal quantum number 
n can be nulled using a (7c(2, 2, ujd) that is positive i.e. by 
setting u!d slightly greater in frequency than the n n-fl 
transition frequency uje,g, so that (5 > 0. 

To null the polarizability difference between two circu¬ 
lar Rydberg states n and n -I- 1, it is possible to choose a 
dressing frequency that increases the Stark shift of the 
lower-n state, while decreasing the Stark shift of the 
upper-n state by working at a frequency intermediate 
between the n -I- 1 — >■ n -I- 2 (lower) and n — )■ n -I- 1 (upper) 
transition frequencies. 

Given the flexibility in the choice of the dressing fre¬ 
quency, how should it be chosen? In Ref. [15], it was 
shown that a dressing frequency could be chosen to ob¬ 
tain dipole nulling and set the differential ac Stark shift 
between the states equal to zero (known as a “magic 
wavelength” in the optical domain [26]). This has the 
clear advantage that spatial inhomogeneities in the dress¬ 
ing field over a sample will have a reduced effect on the 
transition frequency. (The inhomogeneities will however 
lead to an reduction of the efficacy of the nulling.) 

Unfortunately, for the e —> g circular-circular transi¬ 
tion considered here, the (7e,g(0,2) coefficient that deter¬ 
mines the differential ac Stark shift is always positive and 
shows no zero in the frequency range of positive (7e,g(2,2) 
suitable for nulling (see Fig. 2). (The upper e state is 
blue-shifted and the lower g state is red-shifted.) 

However, it is possible to minimize the ac differen¬ 
tial Stark shift: as (7e,g(0,2) scales like 1/(5^, whereas 
(7e,g(2, 2) scales like 1/5^, a lower Stark shift under polar¬ 
izability nulling conditions can be obtained with smaller 
,5. 

As 6 is decreased for a lower ac Stark shift, one must 
consider that sources of the dressing field will not be 
perfect, with spectral impurity of increasing magnitude 
closer to the carrier frequency (see for example Ref. [27]). 
This noise can drive the resonant transition [15] — an 
undesirable situation. However it is also true that to 
obtain nulling, a weaker dressing field amplitude is re¬ 
quired as the dressing frequency approaches the reso¬ 
nance. An optimal compromise between these compet¬ 
ing effects will depend on the spectral noise properties of 
the available dressing source. Somewhat arbitrarily, we 
choose to study (5/27r = —100 MHz, corresponding to an 
offset from the zero dc-field resonance beyond which the 
phase noise density of typical microwave sources is white 
[27].^ 

Figure 3 shows some relevant resonances as a func¬ 
tion of dc field. The choice of 5/2'it = —100 MHz avoids 
multi-photon resonances, such as those that occur at 
« 49.65 GHz. 

With a detuning of 100 MHz to the red of the n = 
50 51 transition {6/2tt = —100MHz), summation 

of the relevant terms of Eq. 11 gives (7e,g(2, 2, We,^ + 
d) « 9.968 Hz/(V/m)"* (application of the simple Eq. 14 
gives (7e,9(2,2, -t- (5) ~ (7^(2,2, -t- <5) ~ 
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FIG. 4. Influence of dressing fields on the dc elec¬ 
tric field dependence of the \n = 50 ,1 = 49, m = 49) ^ 

|n = 51, £ = 50, m = 50) (e g) circular to circular state 
transition with parallel dc and dressing fields. These are 
calculated by diagonalization of the Floquet Hamiltonian for 
varying Fdc- For comparison, the transition frequency with 
no dressing fields is shown. In both the polarizability- and 
dipole-nulled cases the dressing held frequency is 100 MHz 
to the red of the zero ac/dc field eg transition. For polar¬ 
izability nulling, Fac ~ 3.195 V/m; and for dipole nulling, 
Fac ~ 2.718 V/m. The ac Stark shifts in zero dc electric 
field of the polarizability-nulled and dipole-nulled cases are 
~ 0.273 MHz and ~ 0.198 MHz respectively. 


9.90Hz/(V/m)4). With Ce,g{2,Q) « -25.44Hz/(V/m)2 
this predicts that Fac ~ 3.195 V/m will eliminate the 2nd 
order Stark effect — so-called polarizability nulling. 

Polarizability nulling with these parameters can be ver¬ 
ified by diagonalization of the Floquet matrix as a func¬ 
tion of applied dc field i.e. a non-perturbative calculation. 
Figure 4 shows the eg energy difference as a function of 
dc electric field both with and without the dressing field. 
The dipole nulling case shown in the figure will be dis¬ 
cussed shortly. 

It is also possible to null the absolute polarizabil¬ 
ity of either the g or e state (but not both simultane¬ 
ously). More specifically, considering the g state, but 
with a blue detuning of LOdf^n from We.gf^ir of d/2Tr = 
-1-100 MHz, we find that summation of the relevant terms 
of Eq. 11 gives Cg{2,2,ujd) ~ 9.83Hz/(V/m)^. With 
an absolute polarizability corresponding to Cg(2,0) « 
—203.24Hz/(V/m)^, a field strength of Fac ~ 9.09V/m 
is predicted for polarizability nulling. A Floquet diago¬ 
nalization at this Fac shows that quadratic variation of 
the dressed g state energy with dc field is reduced to « 
— 17.3Hz/(V/m)^. Adjusting Fac to 9.54 V/m brings the 
magnitude of this quadratic variation to < 1 Hz/(V/m)^. 

The preceeding discussion demonstrates that it is pos¬ 
sible — in principle — to completely null either the ab¬ 
solute polarizability of a circular Rydberg state, or the 
differential polarizability between two Rydberg states. 
However, we have assumed that the dc and ac fields point 
in the same direction. In practical work with circular 


Rydberg states [1], small uncontrolled electric fields with 
random directions make it necessary to stabilize circular 
Rydberg states against mixing with lower £ and m states 
[28]. For this purpose, a small dc electric field is used to 
remove the energy degeneracy. For example, the exper¬ 
imental work of Ref. [29] used a stabilizing dc bias field 
of 36 V/m. 

To understand the practical aspects relating to non¬ 
parallel fields we will first consider how dipole nulling 
can be obtained in small dc electric fields parallel to the 
dressing field (viewing this as a perturbation from the 
polarizability nulling situation), and then investigate the 
influence of 1 ) fluctuations in electric field transverse to 
the deliberately applied electric field, and 2 ) the conse¬ 
quences of imperfect alignment of the dc bias and dress¬ 
ing fields. 

Dipole nulling at small non-zero dc fields can be un¬ 
derstood using a simple theory. Terms with coefficients 
Cfc(4,0), and Ck{4,,i,ujd) (with i = 2,4,...) in the gen¬ 
eral expansion AEk = '^^jCk{i,j){Fdcy{Fac/2y de¬ 
scribe the quartic variation of energies with dc field. In 
the presence of perfect polarizability nulling, these pre¬ 
dict the range of dc fields in which nulling is effective. 
Consider that a quartic of the form y = oq a 2 X^ + 
a^x'^ has local extrema at a; = 0, ±y/— 02 / 204 . Like¬ 
wise, if polarizability nulling is not perfect {Ce,g{2,£)) + 
Ce,g{2, 2,(jJd){Fac/2)^ y 0 ), the combination of the resid¬ 
ual quadratic variation and quartic terms can lead to a 
local extrema (dipole null) at a non-zero dc electric field. 

For |Ce, 3 ( 4 ,2,u;d)(Eac/2)^| > \Ce,g{^,^)\ (see below) 
the location of the dipole null can be estimated as: 


Fdc 


-Ce,g{2, 0) - Ce,g{2, 2, LOd){Facl2f 
2Cc,g{4.,2,LOd){Fac/2Y 


1/2 


(15) 


Near the We,g resonance, the value of Ce^g{'^,2,0Jd) 
can be estimated using similar considerations as for 
Cc,g{2,2,ujd), again due to the energy degeneracies of 
high-£ states. In fact, the results of Eq. 14 generalize to 
(in atomic units): 


Cc(t 2, OJc^g ^) 


\{uh\F^ kb)r \ {wh\ti^ |cb)P 

Jl-l-l 


for i = 0,2,4, 6 ,.... We can verify that for the n = 
50 ^ 51 system, with 5/2t: = —100 MHz, Fac ~ 
3V/m and Ce,g{‘i^,2,0Jd) « —( 70 ( 4 ,2,We.g + S) that 
|(7e,g(4,2,Wd)(Fic/2)^| > |(7e,g(4,0)| and thus Eq. 15 
is reasonable. To put a dipole null at 40 V/m (ap¬ 
proximately the same dc field used in Ref. [29]) with 
(5/27r = -100 MHz, Eq. 15 predicts Fac = 2.798 V/m. 
(Note that this Fac has a magnitude slightly smaller than 
that required for polarizability nulling since Ce,g{2,0) < 
0, Ce,g{2,2,u}d) > 0 and (7e,g(4,2, Wd) > 0.) A non- 
perturbative Floquet calculation using this value of Fac 
exhibits a null at 36.6 V/m. 

It is straightforward to improve on the estimate of 
Eq. 15 by considering all terms with coefficients of the 
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iorm Ce,g{i,2,ujd)- Differentiating 

Ee,g = CeA‘^,0)Fl+ Y. Ce,g{^,2,^d){FdcnFac/2f 

i=0,2,4,... 

with respect to Fdc, setting this derivative equal to 
zero, and rearranging gives an expression for Fac- With 
Ce,g{i,2,ujd) ~ —Cc{i,2,ujd) as given by Eq. 16, sum¬ 
mation of the resulting arithmetico-geometric sequence 
gives (in atomic units): 

p {-J}3/2{_C^_^(2,0)}1/2 

“ KublM.kb)! K^cblM.lcb)! 

E'dcK'abI ^dz I'M^b)! ^ 

to obtain a null at Fdc- For Fdc = 40 V/m, this equation 
predicts Fac ~ 2.715 V/m. A non-perturbative Floquet 
calculation using this Fac gives a null at Fdc = 40.1 V/m. 
With Fac ~ 2.718 V/m the null is within 10“^ V/m of 
40 V/m — see Fig. 4. 

These perturbation theory results for both dipole and 
polarizability nulling are easy to apply and develop intu¬ 
ition. However, it is prudent to confirm their predictions 
by diagonalization of the Floquet Hamiltonian, partic¬ 
ularly in situations with larger Fdc and Fac, where the 
results of perturbation theory are less likely to be valid. 
For example, the Appendix discusses the dipole nulling 
situation of Ref.’s [11] and [12], where Fdc ~ 400 V/m 
and Fac ~ 88 V/m — both much greater than for the 
dipole null at Fdc = 40 V/m. With 6 « —957 MHz and 
Fac ~ 88 V/m, Eq. 17 predicts a Fdc ~ 265 V/m. Fig¬ 
ure 5b shows a dipole null at 328 V/m, in only moderate 
agreement with Eq. 17. 

The qualitative difference between uje,g for the larger 
FdcS in Fig.’s 4 and 5 (including the presence of the sec¬ 
ond null in Fig. 5) can be explained through examination 
of Fig. 3 where the two dressing frequencies are marked 
with horizontal lines. The upwards trend of the nulled 
cases in Fig. 4 is due to the avoided crossing between 
Iff) ® Iff = 0) and |n = 51, ni = 0, m = 49) (8) |ff = —1), de¬ 
creasing the dressed g energy, which increases the uje,g 
frequency. However, at the significantly lower dressing 
frequency of Fig. 5, there is an avoided crossing between 
Iff) <8 Iff = 0) and |n = 52, m = 2, m = 49) 8 |ff = —2) (a 
two-photon resonance) that pushes the dressed g energy 
to higher energies, decreasing the We,g frequency with in¬ 
creasing Fdc- 


B. Circular Rydberg states with non-collinear dc 
and dressing fields 

It is impossible to perfectly align the dressing and dc 
bias fields directions, and field fluctuations may occur in 
directions other than the bias field. For these reasons it 
is necessary to consider the influence of non-collinear dc 
and dressing fields on dipole nulling. 


As illustrated in Fig. lb, the circular states in zero 
field are energy degenerate with states of the same n but 
lower-m. In the case of cylindrical symmetry (with z the 
axis of symmetry), states of different m are not coupled 
and thus non-degenerate Rayleigh-Schrodinger perturba¬ 
tion theory may be applied, as in the previous section. 
But when the dressing and bias field are no longer orthog¬ 
onal, coupled degenerate states are now involved and the 
non-degenerate approach is no longer applicable. In this 
case we have taken the approach of direct diagonaliza¬ 
tion of the Floquet Hamiltonian. In this situation, it is 
critical that for the basis set parameter 5m > 1. How¬ 
ever, 5m = 1 is entirely adequate for examination of the 
second order effects to be discussed here. 

For concreteness, we consider the ac field to be aligned 
in the z direction and the dc field to deviate from this 
direction. In particular, the scenario of interest is that 
Fdc and Fac are initially perfectly aligned, and that the 
dipole-nulling situation has been achieved i.e. there is 
no first order dependence on deviations of the magni¬ 
tude of the dc field. But in addition to small fluc¬ 
tuations of the field in this direction, we must con¬ 
sider field fluctations in the orthogonal direction: Fdc = 
(Fdc,II + AFdc,i|)z -k AFdc,±x where Fdc,\\ is the dc bias 
field discussed in the previous section. Considering the 
case of Fig. 4, with Fdc^\\ = 40 V/m, we find that to 
2nd order, deviations in the transition frequency are 
given by Aa;e,g/27r « A:j_(AFdc,_L)^ + fc||(AFdc,||)^ (sym¬ 
metry precludes the cross terms) with k\\ and shown 
in the first row of Table I. The magnitude of these 
k's may be compared to the polarizability difference in 
zero dc field with no dressing fields present Ce,g(2,0) « 
—25.44Hz/(V/m)^. The relatively strong dependence of 
LOe^g on the transverse field fluctuations (i.e. large fcj_) 
possibly limits the effectiveness of dipole nulling. For ex¬ 
ample, near a polycrystalline metal surface, patch fields 
show comparable fluctuations in the directions both nor¬ 
mal and orthogonal to the surface [30]. However, excess 
technical noise contributions are often larger in a specific 
direction [31]. 

In the absence of a general theory, we have per¬ 
formed a survey of fey and fc_L for dressing frequencies 
where Ce,g(2, 2, w^) is positive (see Fig. 2). For a given 
LOd we have computed the required ac field strength 
to put dipole nulls at specified FdcH- The results are 
summarized in Table I (which includes the scenario of 
Ref.’s [11, 12]), indicating that there is a slight trade¬ 
off between the variations observed in the two directions. 
Determining whether nulling provides an advantage over 
the normal situation with a bias field — with a linear 
dependence of the transition energy for fluctuations in 
one direction — depends on the specific magnitudes of 
the fluctuating fields one is trying to mitigate against. 

In addition to random field fluctuations, we must also 
consider systematic misalignment between the bias and 
dressing field directions. Small angular displacements of 
these two fields above a threshold value of between 3 
and 4 degrees can completely eliminate the dipole nulls 



at Fdc ~ 328 V/m and w 400 V/m, as shown in Fig. 5 
(curve b) of the Appendix (the inflection point between 
these two extrema vanishes). This sensitivity to an ex¬ 
perimental imperfection is of concern for implementation. 
On the other hand, the dipole null at Fdc ~ 40 V/m is 
Fig. 4 is more robust: the null remains for angular field 
separations up to at least 10 degrees. 


C. Rb Rydberg s-states with non-collinear dc and 
dressing fields 

In light of the high sensitivity of the dipole-nulled 
circular-circular transitions to transverse field fluctua¬ 
tions — as discussed in the previous section — we now 
consider the transverse field sensitivity of the dipole- 
nulled Rb, 49 si/2 —>■ 48si/ 2 two-photon transition stud¬ 
ied previously [15]. 

Without any dressing field applied, the 49si/2 — 48si/2 
transition frequency has a quadratic sensitivity to dc field 
characterized by < 7498 , 485 ( 2 , 0 ) « —294Hz/(V/m)^. Se¬ 
lection of the dressing field frequency and amplitude re¬ 
quired for polarizability nulling is similar to that pre¬ 
sented for the circular case. A plot of < 7495 , 488 (2, 2, w^) 
allows identification of dressing frequency ranges where 
nulling is possible (< 7495 , 485 ( 2 , 2 , w^) > 0 , so that 

<7495,485(2,0) -f <7498,485(2,2, Wd)(Aac/2)^ = 0 may be 
satisfied). Unlike in the circular case, it is possible 
within at least one of these frequency ranges to satisfy 
< 7495 , 485 ( 0 ,2, Wd) = 0. This eliminates the 2nd order dif¬ 
ferential ac Stark shift and makes the transition less sen¬ 
sitive to variations in ac field strength (e.g. due to spatial 
inhomogeneities). 

Plots of < 7495 , 485 ( 0 ,2, Wd) and < 7495 , 485 ( 2 , 2, Wd) vs. Wd 


TABLE I. The sensitivities of the dressed circular-circular 
g —>■ e transition frequencies to variations in the dc field in 
directions both parallel (A:||) and perpendicular (fcu) to the dc 
bias field and dressing field. For a given dressing frequency 
a;d/27r, the dressing field amplitude Fac required to put a 
dipole null at Fdc,\\ has been computed (Fdc,|| and uJdl'2iT are 
labelled as points Fig. 3). The first two lines of this table cor¬ 
respond to situations discussed extensively in the main text 
{S/2n — —100 MHz) and the Appendix (a comparison with 
Ref.’s [11, 12]). Except where noted, the values of cJd/27r and 
Fdc,\\ are exact, whereas the others are displayed to a precision 
consistent with their estimated uncertainty. 


Wd/27r 

Hz 

Fdc, II 

(V/ 

Fac 

m) 

fell 

k± 

Rz/iV/mf 

« 51.00 

40.0 

2.7 

18.3 

5820 

« 50.14 

400.0 

88.0 

-15.9 

6700 

50.35 

200.0 

62.1 

6.8 

16900 

50.35 

500.0 

37.2 

73.1 

969 

49.50 

40.0 

70.4 

5.7 

368000 

49.00 

200.0 

69.7 

5.2 

15000 

49.00 

500.0 

49.3 

48.4 

1170 

48.50 

40.0 

17.7 

1.5 

47000 


are similar to the plots presented in Fig. 1 of Ref. [15] 
for Aa{ujd) and Aj3{ujd), with the correspondence: 
< 7495 , 485 ( 0 ,2, Wd) « —Aa{ujd) and < 7495 , 485 ( 2 , 2 , w^) « 
2A/3{uJd)/Fdc,o where Fdc,o = 100 V/m. In Ref. [15], 
Aj5{u!d) is computed using third-order perturbation the¬ 
ory. The “unperturbed” eigenstates are those obtained 
by diagonalization of the Stark Hamiltonian in a non-zero 
dc electric field {Fdc,o = 100 V/m). Deviations from this 
dc field together with the ac field constitute the pertur¬ 
bation. Likewise Aa{ujd) was computed in a dc field of 
100 V/m and thus is approximately equal but not iden¬ 
tical to C495,485 ( 0 , 2 , CUf/). 

For consistency, we choose a;d/27r = 38.465 GHz 
as in Ref. [15], which is slightly displaced from the 
< 7495 , 485 ( 0 , 2 , Wd) = 0 condition, but gives roughly 

zero differential ac Stark shift at a dc field of « 
lOOV/m. With this Wd, we compute < 7495 , 485 ( 2 , 2, Wd) « 
7.18Hz/(V/m)^, and thus Fac ~ 12.79V/m should lead 
to polarizability nulling (using < 7495 , 485 ( 2 , 0 ) mentioned 
previously). A non-perturbative Floquet calculation with 
a basis set suitable for the Rb states (see discussion in 
Section HB) confirms that the dc polarizability is re¬ 
duced to w 10Hz/(V/m)^ under these conditions. The 
value of Fac can be tweaked to completely eliminate this 
polarizability, or to place a dipole null at a specific field, 
as in the circular-circular case. With Fac ~ 13.01 V/m, 
a dipole-null can be placed at Fdc = 110 V/m. 

The effect of non-parallel Fdc and Fac can be de¬ 
termined using diagonalization of the Floquet Hamilto¬ 
nian. The effect of transverse field fluctations can be 
characterized in a similar way to the circular-circular 
case: for a dipole null at Fdc,|| = 110 V/m, we have 
k\i « 70Hz/(V/m)^, k± w 34Hz/(V/m)^, both of 
which are relatively small compared to < 7495 , 485 ( 2 , 0 ) « 
-294Hz/(V/m)2. 

There is a reduced sensitivity to transverse fluctua¬ 
tions for the dipole-nulled Rb 49si/2 — 48si/2 transition, 
as compared to the circular-circular case (i.e. lower k± for 
the Rb case). This is not surprising given the isotropy 
of the unperturbed s-states compared to circular states. 
We have also found a relatively small sensitivity to trans¬ 
verse fluctuations in the case of a dipole-nulled s—p tran¬ 
sition between triplet Rydberg states of He (manuscript 
in preparation). 

IV. CONCLUDING REMARKS 

Summarizing: 

1. Circular Rydberg states can be polarizability 
nulled, in a similar manner to low-f Rydberg states 
[15]. Both the absolute and the differential po¬ 
larizability between consecutive n, n -I- 1 circular 
Rydberg states can be nulled. The corresponding 
perturbation theory description is relatively simple 
compared to the case of low-f states, where nu¬ 
merical summation of the perturbation series is re¬ 
quired. 
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2. By considering deviations from the polarizability¬ 
nulling situation, the parameters required for the 
dipole-nulling of circular Rydberg states [11, 12] 
can be calculated in a straightforward manner. 

3. Under the conditions of circular state dipole¬ 
nulling, large energy shifts are observed for small 
dc electric field perturbations transverse to dress¬ 
ing and dc bias fields (both parallel). This may 
limit the usefulness of dipole nulling depending on 
the specifics of the electric field fluctuations. 

4. Dipole-nulling of transitions between low-f Ryd¬ 
berg states — such as the 49si/2 — 48si/2 transition 
in Rb — do not show the same high sensitivity to 
transverse fields as circular states. 

Any planned experimental work on the reduction of elec¬ 
tric field sensitivities of circular Rydberg states should 
take the third point into careful consideration. 

Low-ti Rydberg states remain interesting targets for 
the reduction of dc electric field sensitivity using dress¬ 
ing fields — without large transverse field sensitivity. 
Particularly interesting is the possibility of dipole-nulling 
single-photon transitions between Rydberg states of rel¬ 
atively low-^ optically accessible Rydberg states. 
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Appendix A: Comparison with previous dressed 
circular Rydberg atom calculations 

The techniques that we have used in the main body 
agree qualitatively but not quantitatively with previous 
circular Rydberg state dressing calculations [11, 12]. Us¬ 
ing Ref. [32] we have determined the approximations re¬ 
quired for us to obtain agreement with this earlier work. 
As this information is not widely available, we describe 
these approximations and the resulting deviations from 
our approach. 

References [11] and [12] show that a linearly polarized 
non-resonant dressing field can suppress the difference 
between the permanent electric dipole moments of two 
dressed circular states, n = 50, £ = m = 49 ((/ hereafter) 
and n = 51, ^ = m = 50 (e hereafter) in a dc field 
of 400 V/m. In other words, the transition frequency 
between these dressed states is insensitive to first-order 
dc electric field fluctuations (“dipole-nulled”). 

The dipole nulling condition can be obtained using 
different combinations of dressing frequencies and field 
strengths. In Ref.’s [11, 12] the field frequency and am¬ 
plitude are characterized using the near resonant tran¬ 
sition: 5 —>■ f where i is the n = 51, ni = 0, n 2 = 1, 


TO = 49 parabolic state (see Fig. lb). In particular, 
the dressing field amplitude Fac is chosen such that 
Fac\{i\Ig) \/h = 200MHz (where h is Planck’s con¬ 
stant). The transition dipole moment is computed us¬ 
ing the undressed zero dc-field states (i.e. the parabolic 
states are not exact energy eigenstates once a dc field is 
applied — see Section 52 of Ref. [17]). The formulae of 
Ref. [17] give ] {i\g,z \g) \ ~ 177.6 (in atomic units) and 
thus Fac ~ 87.99 V/m. 

The dressing field frequency required for nulling with 
this Fac is specified in Ref.’s [11, 12] as a detuning relative 
to the undressed g ^ i resonance frequency at 400 V/m. 
This resonance frequency is calculated using diagonaliza- 
tion of the Stark Hamiltonian in a finite basis of spher¬ 
ical states \n'£'m'), that includes all valid states with 
m' = n — 1, n' ranging over n, n -|- 1, ..., n + Sn, and £' 
ranging from to^ to n' — 1, where n = 50. The parameter 
Sn controls the size of the basis set and subsequent accu¬ 
racy of the calculation (see below); following Ref. [12] we 
used Sn = 5 (M = 6 in their notation). The resonance 
frequency at 400 V/m is found to be « 50.702299 GHz. 
With a red detuning of 0.555907GHz [12], the actual 
dressing frequency is « 50.146392 GHz [25]. 

This dressing field frequency and the above-mentioned 
amplitude can be used to compute the difference between 
the dressed g and e states as a function of dc field (see 
Fig. 5, curve a). As in the main text, the Floquet method 
is used, with sidebands up to and including ±4 added 
{Sq = 4 in the notation of the main text; iV = 4 in 
the notation of Ref. [12]). As the two circular Rydberg 
states are not coupled by either the dc or ac electric field 

— assuming both fields have the same linear polarization 

— their Hamiltonians can be diagonalized separately. In 
both cases a total of M = 6 different n manifolds are 
used, starting with the lowest n required (n = 50 for g 
and n = 51 for e). 

When calculated in this manner, the energy difference 
between the g and e states does exhibit dipole nulling, 
but at a dc electric field of « 417 V/m, shifted with re¬ 
spect to the value of 400 V/m reported in Ref.’s [11, 12]. 

We find that by either changing the coupling from 
200 MHz to « 194.79 MHz or by changing the dressing 
frequency from « 50.146392 GHz to « 50.142205 GHz 
(see curve b of Fig. 5), the dipole nulling condition can 
be obtained at 400 V/m. 

Thus the procedure of our main text gives results that 
are inconsistent with those in Ref.’s [11, 12] when we 
use the same dressing field parameters. Examination of 
Ref. [32] suggests that this difference is due to the follow¬ 
ing approximation scheme (HM hereafter): ac field cou¬ 
plings are computed in the parabolic basis [17], assuming 
that the zero field parabolic states are the eigenkets of 
the non-zero dc electric field Hamiltonian. (Expanding 
the zero field parabolic states in terms of the spherical 
states [33] allows computation of these couplings using 
Eq.’s (63.2) and (63.5) of [17] — with the correction men¬ 
tioned in Section H.) The energies computed using this 
basis after adding the dc field, but before including the 
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FIG. 5. Transition frequencies between the dressed g and e 
states as a function of dc electric field. The frequencies are 
relative to the zero dc and ac electric field transition frequency 
AE/h where AE — 0.5 (1/50^ — 1/51^) (in atomic units). 
Curve (a) uses our calculation procedure with the dressing 
field parameters of Ref. [12]; curve (b) corresponds to our 
calculation procedure with the same dressing field amplitude 
as Ref. [12], but with a frequency chosen to place a null at 
400 V/m; curve (c) is a calculation using the HM procedure 
(see text) with the dressing field parameters of Ref. [12]. The 
inset compares (b) and (c) around the dipole nulled condition, 
relative the respective dressed transition frequencies at Edc = 
400 V/m. 


ac field, are assumed to be the 2nd order perturbation 
energies (Eq. 52.3 of Ref. [17]), except for the g, i and 
e states; in this case the energies are obtained from a 
full diagonalization in a M = 6 sized basis set (using the 
notation of Ref. [12]). The essence of this approximation 
is that the ac field couplings do not consider that the dc 


field can mix states of different n. 

With the HM approximation scheme and the parame¬ 
ters for the dressing field amplitude and frequency com¬ 
puted from Ref. [12] (see above and [25]) we have ver¬ 
ified that dipole nulling occurs at a dc electric field of 
400.024V/m, consistent with Ref.’s [11, 12] (see Fig. 5 
curve c). As expected, the HM approximation scheme 
agrees with the more complete calculation as the dc field 
goes to zero (then the ac couplings become exact). 

Figure lb of Hyafil et al. [11] contains a plot with a sim¬ 
ilar axis to the inset of our Fig. 5. When their Fig. lb 
is digitized we find agreement with the HM approxima¬ 
tion scheme to within 0.03 V/m and 0.25 Hz from 399 to 
401 V/m. 

Mozley et al. [12] discuss the size of basis set required 
to reach a particular accuracy. They found that for Sn > 
5 and Sq > 4, a. subsequent increment in either Sn or 
Sq did not shift the dressed energies by more than 1 Hz. 
Because it might be possible that the HM approximation 
scheme significantly influences the size of the basis set 
required, we have also checked the values of Sn and Sq. 

The Sn required depends on the dc electric field. In 
zero dc field we find that i5n = 5 is adequate for both our 
and the HM method; but at 400 V/m it is necessary to 
go to Sn = 7 to meet the 1 Hz criteria for both methods. 
However, this 1 Hz convergence criteria is less stringent 
when applied to the energy differences as compared to 
the energies, which tend to shift similarly as the basis set 
is increased. If we consider the dc field at which nulling 
is obtained, we find that curve b of the inset of Fig. 5 is 
shifted by only 0.001 V/m between Sn = 5 and Sn = 7. 

In summary, we conclude that good quantitative agree¬ 
ment is obtained with the results of Ref.’s [11, 12] when 
we use the approximations described in Ref. [32] (the HM 
scheme). This comparison has been presented in the in¬ 
terest of reproducibility; we do not advocate use of the 
HM scheme as we have found it more difficult to imple¬ 
ment and of no advantage in computational efficiency. 
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